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Abstract. We investigate the radiative decays of the </> meson to the scalar mesons ao(980) and /o(980). 
We demonstrate that, contrary to earlier claims, these decays should be of the same order of magnitude 
' for a molecular state and for a compact state and, therefore, the available experimental information is 

, consistent with both, a molecular as well as a compact structure of the scalars. Thus, the radiative decays 

. of the (j) meson into scalars establish a sizable KK component of the scalar mesons, but do not allow to 

' discriminate between molecules and compact states. 

^ ! 
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Qh! 1 Introduction 

Oh: 

(D . 

'T^ . It has been claimed for many years that studies of radiative decays (f> 709 (980) — > 77r°?7 and (f> 7/0 (980) '^tt'^tt'^ 
' are a powerful tool to discriminate between various models for the low-lying scalar mesons. The extraction of the (jj^ao 
and (/)jfo coupling constants from the data is not a straightforward task (see ^), but it is a common belief that, with 
data accurate enough, radiative decays would reveal the nature of the lightest scalars. 

5^ , The simplest mechanism for these radiative decays assumes that the ao and /o are "^Po quarkonia, and the decays 
proceed via a quark loop. Nevertheless, with the (/)-meson being mostly an ss state, this mechanism cannot be re- 
sponsible for the decay 4> 700, since, in the quarkonium picture, the oq is an isovector state made of light quarks. 
Similarly, only /o(ss) can be produced via the quark loop mechanism and, if so, the subsequent decay /o tt^tt^ is 
suppressed by the OZI rule. On the other hand, as both /o and uq are close to the KK threshold and are known to 
couple strongly to this channel, one expects that the radiative decay mechanism via charged kaon loop should play 
an essential role, as it was suggested in Refs. The existing data on cf) radiative decays [^||^][7j support this 

expectation, as is shown in detail in Ref. 

The latter observation does not mean per se that the quarkonium assignment for ag and /o is excluded by the data. 
It only means that the strong coupling to the KK channel, together with the threshold enhancement phenomenon, 
makes the kaon loop mechanism dominant. However, the strong coupling to KK implies that the KK component in 
the wave functions of these mesons should be large, and recent studies 9 based on the analysis of near-threshold data 
confirm this. A large KK admixture should be reflected somehow in the radiative decay amplitude. 

In Ref. PI it is claimed that there should be a strong suppression of the <j> 7/0 / ao branching ratio for the scalars 
in case they are loosely bound molecules as compared to pointlike scalars that correspond to compact quark states, 
(10~^ vs 10"^). A study by Achasov et al. 0], where the finite width of scalars was taken into account, arrived at the 
same conclusion. Thus, the authors of P] and ^ stress that data for this branching ratio should allow to prove or 
rule out the molecular model of the scalars. However, no such suppression was found in recent kaon loop calculations, 
Refs. |10II11II12| . where the scalars were considered as dynamically generated states, i.e., as molecules. The aim of 
the present paper is to demonstrate explicitly the implications of a molecular structure of scalars on the radiative (j) 
decay. In the course of this we can demonstrate what went wrong in the analysis of Ref. [3] and confirm the results of 

Refs. uniiniin]. 
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2 Point-like scalars 

To simplify the situation we work with stable scalars — the generalization to a more realistic case is straight forward 
and should not change the conclusions; we comment on what is necessary for this generalization in what follows. The 
current describing the radiative transition between the vector meson and a scalar meson S, in the kaon loop model, 
is written as |13II14| (see [3] for notations) 



9<p9s 



■I{a,b)[e^{p ■ q) - p^{q ■ e)], 



(1) 



where p and q are the momenta of the 0-meson and the photon, respectively, rriK is the kaon mass, and gs are the 

(jjK^K^ and SK^K^ coupling constants, Si, is the polarization four- vector of the 0-meson, a = — ^, and b — — #- (in 

case of an unstable particle produced rrig is to be replaced by the invariant mass squared of the decay products). The 
amplitude |^ is transverse, M^q^, = 0, and is proportional to the photon momentum. 

For the pointlike model of the scalar mesons the function I {a, b) was calculated in Refs. j2E|- It is given by 



I{a,b) = 




5(a) 



V4a- larcsin(2^). 



In(^) 



a > 
a < 



9 



(2) 



7/± = ^ (1 ± VT^) ■ 

Note that the integral /(a, b) remains finite in the limit a ^ b. 

To arrive at the formula consider the sum of the graphs depicted in Fig. l(a)-(c), where the appearance of the 
graph 1(c) is a consequence of gauge invariance, since the (j) KK vertex is momentum-dependent. The current in 
Eq. (P) is given by = eg^gse^Jf^^,, with 



where 



d*k 



J — ji"-) 4- jib) I 7(c) _ O j(a) I 7-(c) 

{2k-p)^{2k~q), 



(27r)4 [fc2 



m? + iO] [{k — q)"^ — + iO] [{k — p)^ — + iO] ' 
d*k 1 

(27r)4 [fc2 - m2 + iO][{q + k - pf - m? + iQ] ' 



(3) 

(4) 
(5) 



and m = mx- 

Since gauge invariance demands the structure of the integral ijjjl to be 



Jt,y = J[p,^qt, - [p ■ q)9t^A 



(6) 
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the strategy applied in Ref. jSj is to read off the coefficient of the p,yq^ term, coming entirely from the integral 
and to restore then the coefficient of the g^,y term with the help of Eq. This allows the authors to deal with a 
finite integral and thus to bypass the problem of treating the divergent parts of the loop integrals l@J, (jS)). However, 
as we shall see below, the divergent pieces cancel and the sum of diagrams given in Eq. Q is finite jl5| . 



To see this we decompose the expression for jjn) as 



where 



1 /"^ . L l-az(l-z)^ l-bz{l-z) 



2n^m? \(a — b) 



dz 







1 - z - 
dz 



z(a — b) 



1 



z 



z)a 



Here and in what follows we consider the case of > 2m, ms < 2m. In addition 

"2 



az(l — z) 



(7) 
(8) 



j: = 



167r2 



- 7B - In 2 



8^ 



dz(l-z)ln[l-6z(l-z)], 



(9) 



where /i^ is the auxiliary mass parameter, the number of dimensions D is equal to 4 — e, and "fE is the Euler constant. 
Similarly, the contact term jSj) can be presented as — 2(;^t^ with 



J' 



167r2 



IE 

e 



In- 



167r2 



dzln[l - - z)], 



(10) 



and, since 



dz{l - 2z) ln[l - bz{l - z)] = 0, 



the structure © is restored. We conclude therefore that, with the proper regularization, the total matrix element is 
finite. It means that the range of convergence of the integrals involved is defined only by the kinematics of the problem. 
In particular, if both masses of the vector and scalar mesons are close to the KK threshold, the integrals converge at 
fco ~ m and for nonrelativistic values of the three-dimensional loop momentum k, |k;| <C m. The nonrelativistic limit 
of the integral I{a, b) takes the form 



/jvi?(a, b) 



+ i- 



7r(.T^ + Sxy"^) 

24(a;2-|-y2)2 ' ■^I2(a;2+y2) 



2^2 ' 



(11) 



where 



y = V(«/4) - 1 x = Vl - (V4) , x,y<$:l. 
Note that, although the expression Hll() contains the factor ^^^^^ ~ T[~b^ does not mean that lNR{a,b) blows up 
in the limit of zero photon energy, w — s- 0. Indeed, the formula Hll() is valid for the scalar meson lying below the KK 
threshold, so one cannot put uj ^ here. If the scalar appears above the kaon threshold, Eq. (jllll is replaced by 

ni {2y + x) 



24 (y + 5)2^ 

so that /(a, b) remains finite in the limit w — s- 0. 



(12) 



3 Introducing the scalar wave function 



When treating the scalar meson as an extended (non-pointlike) object it is not sufficient to insert the corresponding 
form factor into the K^K~S vertex (see 0), but gauge invariance calls for a correction term induced by this additional 
flow of charge. Since only soft photons are involved the needed correction term can be expressed as the derivative of 
the formfactor inserted. Thus we get for the induced vertex 



r,{K+K-Sj) = -2(p+ 



ar(p2,TO2) 

Pu) 



9p2 



(13) 
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where r{p'^,p^_) = r{p^,p'^) parameterizes the momentum dependence of the K^K S vertex, with r{m'^, rn^) — 1. 
Here p+ and are the and four-momenta, respectively. The corresponding extra diagram is depicted at 
Fig 1(d). 

Before proceeding further we note that inclusion of the extra contact vertex (I13|l is a way to insert an ultraviolet 
cutoff in a gauge invariant way. As demonstrated above, the integrals of interest converge already for nonrelativistic 
momenta even for a pointlike vertex, thus it is justified to use nonrelativistic kinematics also when the vertex function 
r is included, as it was done in — one only needs the mild assumption that F decreases faster than 1/k for 
increasing values of its arguments. Then only the positive-energy parts of the kaon propagators are retained, the kaon 
energies are replaced by m, and and 7715 are replaced by 2m, wherever possible. As to the vertex function, in the 
nonrelativistic description the virtuality of kaons is measured by the relative momentum of kaons in the intermediate 
state, so that in the center-of-mass frame of the vector meson [p ~ 0) the vertex function is a function of the 
three-momentum of the outgoing kaons only and thus the spatial loop integrals read 

J,^ = 2ji") + + j[t^ = ~5.k^{a - b)I{a, b; T) + ... , (14) 

when evaluated in the rest frame of the vector meson. Terms that do not contribute to the process of interest are not 
shown explicitly. Note, gauge invariance is enshured by the appearance of the term (a — b) that vanishes for vanishing 
outgoing photon energy. The individual integrals are 



tPk fc,fcfcr(|k-q/2|) 



- m^j (2.)3[^^_^^,o][^^_(k^^^o]^ 
r(c)_ f d^k r{k) 



j{d) 'i_ f d^k kjkk 1 dr(k) 

" 2m^J {27rf Ev-^+iOk dk ' ^ ' 

We assume Ey = my — 2r7i > 0, Es = ms — 2m < for looking at only one kinematic regime is sufficient to make 
our point clear. For more realistic calculations that include the finite width of the scalar mesons we recommend Refs. 
|l()lllllll2j . Performing integration by parts in the integral J^^'^ , one has 

•^.fe - 2m2 j (2^)3 Ev^^+iQ 3m3 j {2nf [Ey - ^ + lOY ' ^ ' 

This trick was used both in Refs. |3] and |llj . 

Let us now assume that F decreases with the range (5 that satisfies the conditions 

(3^->mEv, l3^:>m\Es\. (17) 

With the help of the representation (|16|) one immediately sees that, in the limit [3 00, the divergent terms in J^j,, 
Eq. H14|) . cancel each other and, in the leading nonrelativistic approximation, Ey <C m, \Es\ ^ m, the total matrix 
element does not depend on (3: 

I{a,b;r)^lNR{a,b)- (18) 
We stress that the result (|18|) follows from the nonrelativistic formula H14() . and the only condition needed is H17|) . 
We have repeated the calculation of /(a, 5; F) presented in Ref. jllji with the model form factor /^(k) = /3^/(k +0^). 
The results are depicted at Fig. 2 together with the results of the full pointlike theory. One can see that, in the soft- 
photon limit, there is no considerable suppression of the matrix element due to finite values of (3, down to [3 ^ 0.3 
GeV. The reason for this was discussed above — the integral of Eq. ^ converges for nonrelativistic values of |k|, in 
the soft-photon limit. 

Now we specify the form factor in the molecular model for the scalar mesons. To this end we use the well-known 
quantum-mechanical expressions which relate the KKS vertex and the wave function of the molecule. In the vicinity 
of a bound state the nonrelativistic f-matrix i(k, k', E) takes the form 

t(k,k',i?) = ^MLl, ^(k)=«c/,(k), (19) 

E + e — lO 

where <^(k) is the bound-state wave function in the momentum space, normalized to unity, e = —Es is the binding 
energy, and the Schrodinger equation for the bound state is written symbolically as 

k^ 

— (/.(k) + w0(k) = -e0(k). (20) 
m 
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Fig. 2. The real (the left plot) and the imaginary (the right plot) parts of the function /i = (a — b)I{a, b; F) for P — 0.2 GeV 
(dash-dotted line), /? = 0.4 GeV (dotted line), /? = 0.6 GeV (dashed line), and /? = 0.8 GeV (thin solid line). The result of the 
full pointlike theory is given by the thick solid line. 

The relativistic vertex differs from the nonrelativistic vertex 7 by a kinematical factor (see, e.g., jl6j). 



gsr{k) = (27r)3/V8m2ms 7(k) 



(21) 



where the effective coupling gs is introduced to ensure the normahzation condition r{0) = 1. Using the bound-state 
equation (|20|l . one has, finally. 



gsr{k) = (2^)3/2 



ms\ —+s] 4>{k). 
m 



(22) 



Thus we find that the momentum dependent factor that appears in Eq. (|22|l exactly compensates for the two kaon 
propagator in Eq. ((T3|l . The wavefunction then supplies exactly that piece due to its demanded asymptotics. 

A real molecule is a loosely bound state with a large mean distance between the constituents — much larger than 
the range of the binding force ro . In this deuteronlike case one has 



'me. 



(23) 



Correspondingly, the vertex H22(l does not depend on k, and one can safely use the formulae Q, of the pointlike 
theory with 

(2^)3/2 



gs 



n 

I 9s I — 

-y/imsn, — ~ 32my/me. 
47r 



(24) 



The nonrelativistic expansion of the integral I{a, b) can be used as well. 

So we conclude that the range /3 of the form factor should be identified with the inverse range of the force, f3 ^ l/^o, 
and, if the inequality 

K 

(25) 



Kro ~ < 1 



holds true, the results of the pointlike theory for the radiative 4> decay are valid for molecular model of the 

scalar. In particular, there is no special suppression of the matrix element due to a finite value of /?. 

The latter statement is based on the validity of the inequality (|25|) . What values of (3 would one expect in realistic 
models of the KK molecule? In the meson-exchange models like jl7| it is argued that a strong t-channel force is 
responsible for the formation of scalars. In such a case it is reasonable to identify /3 with the mass of the lightest meson 
exchanged. As there is no pion exchange in the scalar sector, the lightest meson should be the p, which gives for f3 
the value of about 0.8 GeV. In the phenomenological picture of Ref. ^Jj, (3 is taken to be 0.5 0.7 GeV. In the quark 
language, /? is defined by the scale of the internal size of the quark wave function, which also leads to the estimate for 
(3 to be of the order of a few hundred MeV. With such estimates, the inequality H25(l is safely valid for the masses of 
the scalar about 970 ~ 980 MeV. 

The formula (|24|) implies that the vertex gs depends on the binding energy and its value decreases with decreasing 
binding energy. This in turn causes a suppression of the branching ratio when the binding energy tends to zero, cf. 
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Fig. 3. The dependence of the branching Br{(j> "/S) on the mass of the scalar meson in the molecular model. 

Fig. 3. However, for binding energies of typical order of magnitude, for example, e = 10 MeV, Eq. (|24|l yields a coupling 
constant gs of 

^ = 1.12 GeV'^. (26) 
47r 

That corresponds to a branching ratio of i3r((/) jS) « 2.6 x 10"'* which means that there is practically no suppression. 

Nevertheless, we should emphasize in this context that a reliable quantitative calculation of the width certainly 
requires a more realistic approach where it is taken into account that the scalar mesons have finite widths due to 
the presence of the light pseudoscalar channels, and that the quantities that are really measured are the transitions 
(f) — » 77r7r or ^ ^ jTrrj. The impact of finite width effects have been thoroughly investigated by J. A. Oiler ^md also 
by Achasov and Gubin 18 and we refer to their work for details. Here we only want to make the reader aware of the 
fact that due to the proximity of the •jS threshold to the mass of the (/)(1020) resonance even small variations in the 
nominal resonance masses of the scalar mesons have a drastic effect on the available phase space and in turn on the 
obtained results - as it can be imagined from Fig. 3 - unless the finite width of the (/o(980) or ao(980)) scalar mesons 
is considered [1^1115) . 

To take into account finite width effects one has to use the two-channel version of Eq. H19|) from the very beginning 
so that the vertex which appears in the loop integral is accompanied by the vertex that appears in the resonance 
decay, as it is required by the two-channel unitarity condition. If the characteristic scale /3 in this full i-matrix is not 
too small, then the feature that there is no specific suppression due to the molecular structure of the scalar mesons 
will be preserved, cf. appendix A. 



4 Comparison to older work and conclusions 

Our findings are in contradiction with the results of Ref. [H]. The specific model for the KK molecule used there 
was taken from Ref. 19 , which, in turn, is a modification of the approach developed in Ref. |20| and based on the 
quark-exchange picture. The KK interaction employed in Ref. ^19, , was approximated by a local potential of the form 



V{r) — —Vq exp 



1 / r 



ro 



(27) 



with rp = 0.57 fm. This interaction gives e = 10 MeV, so that nro ^ 0.2, and the molecule is rather deuteronlike. 
The wave function was parameterized as 



V'(r) = ( ^ 



1/2 



-fir 



(k) = 



(2/i) 



3/2 



TT (k" + /i2) 



2\2 ' 



(28) 



with fi — 0.144 GeV. This wave function yields a good approximation for the exact wave function, in the momentum 
space (see 0). On the other hand, the wave function H23|l with e = 10 MeV looks very similar, see Fig. 4. 
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Fig. 4. The wave function of the KK system, in momentum space. The approximate solution Eq. (1281 — the curve 1, and the 
deuteronlike wave function, Eq. 1231 — the curve 2. 



So what is wrong with Ref. , and where does the suppression of the radiative decay amphtude come from? The 
answer is rather simple. In Ref. the calculations of the loop integrals were performed with using the wave function 

0(k) = m)^J^ — , (29) 

as a form factor, instead of the correct formula H22|l for the form factor. This led to the result of 4 x 10^'"' for the 
branching ratio (or r{(l) -fS) = 1.7 x 10^'' MeV). The same incorrect choice for the form factor was made in 4 . As 
fj, is as small as 0.144 GeV, no surprise that the suppression found was huge! 

The radiative decay width calculated with the parameterization of the wavefunction (|28|l and the correct formula 
is r{(j) 75) — 2.4 X lO"'^ MeV. It is somewhat large as compared to the experimental result. We want to point 
out, however, that this is primarily due to the not very accurate parameterization. Indeed, the approximation (|28|l 
is definitely wrong for distances beyond the range of the forces, r 3> tq, where the wave function should behave as 
. On the other hand, the deuteronlike wave function is wrong at short distances. It is clear, however, that 



possible contributions to the integral (|14|l coming from short distances correspond to large values of |k| where the 
integrand is suppressed. The value of 1.1 x 10"'^ MeV for the width, obtained in the pointlike theory with a value 
of gs given by Eq. I|2t)|) , is therefore a good approximation to the corresponding width calculated within a molecular 
model [B] . 

In conclusion, there is no considerable suppression of the (p — s- jS width in the molecular model for the scalar 
mesons. As soon as the form factors of an extended scalar meson are treated properly, the corresponding results 
become very similar to those for a pointlike scalar meson (quarkonium) , provided reasonable values are chosen for the 
range of the interaction. We confirm the range of order of lO^'^ 10^** for the branching ratio obtained in Refs. ^1 

unin. 
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5 Appendix: Inclusion of a finite width 

In this appendix we discuss the effect of a finite width of the scalar mesons, due to their decay into two pseudoscalars 
(P1P2), on the total width r{(j) ^ jS). 

The P1P2 invariant mass distribution has the form 



\{a-h)l{aMMKK^p.P.{^P.P.)\\ , 6 - (^^) , (30) 

\ m / V TO / 



amp^p^ 3[2n)°m^ V to / 
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2.56 


3.07 


2.80 








1.22 


1.57 



Table 1. The branching ratio Br((f> — > jS) x 10^; ki,2 are given in MeV. 



where lu = * 2^^^ photon energy and mp^p^ is the invariant mass of the outgoing pseudoscalars. Here the 
range of the force is assumed to be small enough so that one can take the integral /(a, b) for the point-like case, cf. 
Eq. (j23l. 

To account for the finite width of the scalar mesons one is to use the two-channel i-matrix. For the deuteronlike 
case, the amplitude in the KK channel can be written in the scattering length approximation with a complex scattering 
length a^^, ^ 

aKK = — > '«2 > 0, (31) 

Hi + IK2 

for energies around the KK threshold (and energies sufhciently far away from the P1P2 threshold). Then the KK 
P1P2 transition amplitude A squared can be found as 

\AKK^ppAmp,p,)\^ ^ ; ^ = , (32) 

with E = mp-^p^ — 2m. 

In the limit K2 ^ there is no coupling to the P1P2 channel and, for ki > 0, there is a bound state in the KK 
channel with the binding energy e — K\/m. One can readily obtain the total radiative width in this case, which is 
given by the standard formula, 

r((/) -> 75") = 4 2 |(a - b)I{a, b)\ , a=[ 1 , 6=1 1 , oj = ^ , (33) 

with ms = 2m — e and gs defined by Eq. H24() . 

In order to estimate the effect of a finite inelasticity K2, we have calculated the contribution to the total width, 

Ptot ^ [ dmp^P^—^ — (34) 
J dmp^p^ 

from the distribution 1)30(1 integrated over the near-threshold region, 900 MeV < M < m^. The results for the branching 
ratios are listed in Table 1. One can see that the branching ratio remains in the order of 10~^ even for ki = 0, if the 
scale of K2 is around 50 100 MeV. We conclude, therefore, that the results presented in this paper are robust against 
the inclusion of the finite width of the scalar. 

We would like to note here that the above-mentioned scale for K2 is quite natural. For example, as it was shown in 
Ref. (21], the data on the tttt scattering near the KK threshold are, indeed, nicely described in the scattering length 
approximation, with K2 lying in this range (and the ratio ki/k2 being of order unity). 
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